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We demonstrate that migration away from self-produced chemicals (chemorepulsion) generates 
a generic route to clustering and pattern formation among self-propelled colloids. The clustering 
instability can be caused either by anisotropic chemical production, or by a delayed orientational 
response to changes of the chemical environment. In each case, chemorepulsion creates clusters 
of a self-limiting area which grows linearly with self-propulsion speed. This agrees with recent 
observations of dynamic clusters in Janus colloids (albeit not yet known to be chemorepulsive). 

More generally, our results could inform design principles for the self-assembly of chemorepulsive 
synthetic swimmers and/or bacteria into nonequilibrium patterns. 


PACS numbers: 05.65.+b,87.17.Jj,82.70Dd,82.40.Ck 

Active systems, such as suspensions of autophoretic 
colloidal swimmers, motile bacteria, or other self- 
propelled particles, are far from equilibrium even in 
steady state due to their continuous energy expendi¬ 
ture [1, 2]. Unlike isothermal Brownian colloids, motile 
particles can accumulate in regions where they move 
more slowly. They also slow down at high density, creat¬ 
ing a positive feedback loop that can lead ultimately to 
motility-induced phase separation (MIPS) [3, 4]. 

Experiments with artificial self-propelled colloids have 
reported self-organised dynamic clustering, sometimes at 
densities well below that expected to trigger MIPS [5, 6]. 
These “living clusters” seem to reach a limiting size and 
do not coarsen indefinitely: they show microphase sep¬ 
aration rather than macrophase separation as in MIPS. 
So far, the underlying mechanism remains unclear. 

In such experiments, motility is autophoretic: a chem¬ 
ical reaction is catalysed on part of the colloid surface, 
creating a gradient of reagent and/or product that drives 
the particle forward by diffusiophoresis or a similar mech¬ 
anism [7-9]. The same gradients can then cause chem¬ 
ically mediated long-range interactions between the col¬ 
loids, and can also cause rotational torques that bias 
the swimming direction up or down the chemical gra¬ 
dient (an effect known as chemotaxis) [10, II]. This 
fact has suggested a parallel between the experiments 
in [5] and the Keller-Segel (KS) model [12-14] for mi¬ 
croorganisms interacting via chemical signalling. This 
mapping, which assumed that active colloids swim up 
chemical gradients (the “chemoattractive” case), can ex¬ 
plain clustering, but leads to complete phase separa¬ 
tion, rather than a self-limiting cluster size. A combi¬ 
nation of a passive drift up the chemical gradient and 
self-propulsion down it (“chemorepulsion”) might lead to 
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finite size clusters [16]; however, at variance with exper¬ 
iments [5, 15], these should shrink as the self-propulsion 
speed increases [16]. A more general study of chemore- 
sponsive active colloids in the limit of fast chemical dy¬ 
namics suggests a far larger potential for pattern forma¬ 
tion than is predicted by the KS model [10]. 

Here we propose a theoretical framework for active col¬ 
loids that describes at a continuum level not only col¬ 
loidal and chemical densities, but also the local mean 
orientation (“polarization”) of the active particles. We 
thereby uncover new instability mechanisms which show 
that chemorepulsion can lead to clustering and pattern 
formation. These mechanisms were largely ignored in the 
literature so far, perhaps because chemoattraction offers 
a more obvious route to (bulk) phase separation. 

We call our two mechanisms the “Janus instability” 
and the “delay-induced instability”. The physics under¬ 
lying these is indicated in Eigs. Ib and c respectively. Eor 
the Janus instability, active particles reorient and move 
towards a fluctuation-induced minimum in the chemi¬ 
cal density, thereby forming an inward-pointing aster 
(Eig. Ib). Due to local anisotropy of the chemical produc¬ 
tion (inevitable for Janus colloids half-coated with cat¬ 
alyst), this cluster produces a shell of chemorepellent, 
keeping the particles within it together and driving oth¬ 
ers away, so that it cannot grow beyond a certain size. 
In the delay-induced instability, which appears in a dif¬ 
ferent parameter regime, colloids respond differently to 
the same fluctuation. Initially they accumulate at the 
chemical minimum, producing more chemical and lifting 
the concentration there (Eig. Ic). If they reorient slowly 
on the timescale needed to cancel the initial fluctuation, 
the production can overshoot and the minimum becomes 
a maximum before particles finally drift away. This can 
trigger a cyclical instability towards traveling wave pat¬ 
terns. When both instability mechanisms act together, 
we predict “blinking” patterns where clusters continually 
exchange particles and never reach a steady state. 

These chemorepulsive instabilities may shed light on 
the dynamic clustering or microphase separation recently 
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FIG. 1. Schematics: a.) Colloidal Janus particles (spheres) half-coated with a catalytic material (yellow) that produces a 
chemical species (blue); the colloids self-propel by autophoresis in the resulting gradient. Chemorepellent colloids propel from 
high to low concentration of the chemical, chemoattractant colloids vice-versa, b.) Anisotropy-induced Janus instability, which 
can lead to aster-like clusters of self-limiting size, c.) Delay-induced instability, which can lead to travelling wave patterns. 


observed in experiments [5, 6 , 15]. Both instability mech¬ 
anisms induce clusters of a self-limiting area that grows 
with self-propulsion speed, which qualitatively agrees 
with experiments [5, 15]. More generally, chemorepul- 
sive instabilities suggest new design principles for the ac¬ 
tive self-assembly of colloids into spatiotemporal patterns 
with tunable properties. Our key message that chemore- 
pulsion can generate instabilities of uniform states could 
also be relevant for biophysics, where the chemoattractive 
KS instability has long been assumed to drive structure 
formation among microorganisms [17, 18, 20]. 

We describe active colloids (living or artificial) at a 
coarse-grained level, through their density and polariza¬ 
tion fields, p(x, t) and p(x, t). The latter is a local aver¬ 
age of the unit vector describing the propulsive direction; 
this rotates in responds to gradients of a chemical density 
field c(x, t). The colloids self-propel at constant speed 
vq and also have isotropic diffusivity Dp. To represent 
autophoretic colloids (or signaling bacteria), we assume 
that the chemical species is produced by the colloids at 
local rate k^p - with an important, p-dependent correc¬ 
tion addressed below - and decays at rate kd- Therefore, 
we can describe the system dynamics by 

P = -V • (pfop) + -DpV^p 

P = -7P + -DpV^p +Wc-72|P^|P 

c = DcV^c + kop - kdc + kaV ■ (pp) (1) 

Here (3 measures the chemotactic coupling strength; 
when positive, this represents chemoattraction (Fig. la), 
for instance bacteria swimming up food or aspartate gra¬ 
dients [18, 19]. Here however we focus on negative /3, de¬ 
scribing chemorepulsion, as arises for at least some types 
of colloid [10], or for cells fleeing from toxins [18]. In (1) 
the polarization decays locally at a relaxation rate 7 , set 
by rotational diffusion; it also has translational diffusiv¬ 
ity Dp rsj Dp ^ Vq/t [4], which smears out details of p 
on scales below the “run length” vq/^. The term in 72 
describes saturation in p at strong alignment [21]. Fi¬ 


nally, kaV • (pp), where ka has the dimensions of speed, 
describes an anisotropic correction to the isotropic chem¬ 
ical production term {kop), arising whenever the chemical 
is produced by the colloid asymmetrically [22, 23]. 

Rewriting (1) for dimensionless quantities t = kdt, 
X = x^kd/Dp, and setting p = pko\P\vo / {k^Dp)= 
vo/{\/kdDp)p',c = vo\(i\/(kdDp) we show (in SI) that the 
parameter space is spanned by five dimensionless vari- 
ables r = 'y/kd', Vp = Dp/Dp] Vc = Dc/Dp] n = 
kakd/{koVo); Tg = ^ 2 Dplvl. Also s = sign(^) distin¬ 
guishing positive (s = 1) and negative (s = — 1) chemo- 
taxis. In the following, we omit tildes. 

We have solved (1) numerically on a square box of 
side L by finite difference methods. Uniess otherwise 
stated, resuits are obtained for periodic boundary con¬ 
ditions, with as initiai condition a smaii p-perturbation 
of the uniform state (p, p, c) = {po, 0, po)-With chemoat¬ 
traction {s = 1 ), the initiai uniform state is stabie for 
smaii po, whereas for stronger coupiing it is unstabie to 
the formation of dense coiioidai ciusters that co-iocaiise 
with maxima in c (Fig. 2a,b). These dropiets coarsen 
continuousiy to yieid compiete phase separation at iate 
times, aibeit featuring orientationai order in the form of 
a macro-aster (Fig. 2d, inset). This chemicaiiy-induced 
phase separation is weii understood [14]: coiioids swim 
towards high chemicai concentration forming a ciuster, 
which increases chemicai production iocaiiy, recruiting 
further particies, etc.. 

Aithough this feedback ioop is absent for chemorepui- 
sion {s = — 1 ), upon tuning po beyond a certain thresh- 
oid, we observe, strikingiy, that the initiai fluctuations 
ampiify aiso for chemorepuisive coiioids. In marked con¬ 
trast with the chemoattractive case, the resuiting dense 
coiioidai ciusters do not coarsen beyond a characteristic 
size, instead yieiding microphase separation (Fig. 2e-h 
and Supporting Information (SI), Videos 1,4,5). The 
behavior en route to this steady state is compiex; it 
can feature ampiitude osciiiations, or “biinking” (see SI, 
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FIG. 2. Time evolution of the density fields p,c, |p| (insets: p/|p|) for a weakly perturbed uniform initial state (color shows 
value of the fields), (a-d): Clustering and phase separation for chemoattractive colloids (s = l,po = 8). (e-h): Arrested phase 
separation and stationary density pattern for chemorepulsive colloids resulting from the Janus instability (s = — l,po = 10). 
Other parameters: Pp = Pc = = F = 1, r 2 = 10 and Xu = 5; = 1. (Time and space units of r = t/tu and x = xjxu). 


Videos 1,4,5). Blinking clusters dynamically exchange 
particles, before settling down into a stationary arrange¬ 
ment, which normally consists of a hexagonal lattice of 
droplets and an inverted pattern in the chemical density 
(Fig. 2 g). Deviations from the ideal hexagonal structure 
(visible close to the boundaries in Fig. 2) can be more 
or less pronounced, depending on the specific parameter 
choices. As can be seen in Fig. 2 f (inset), chemorepulsive 
colloids point towards the chemical density minima be¬ 
tween the colloidal clusters, and the overall orientational 
pattern consists of stable asters and anti-asters. Param¬ 
eter fine tuning, or choosing no-flux boundary conditions 
instead of periodic ones, can lead instead to permanent 
blinking (SI, Videos 4,5). 

Intriguingly, we also find colloidal waves and travelling 
oscillatory patterns. In particular, upon further increas¬ 
ing the overall density po we observe an amorphous pat¬ 
tern which evolves towards a more regular state, where 
clusters continuously merge, split and decay (Fig. 3 and 
SI, Video 2,3). Eventually, the colloidal and chemical 
density fields may approach regular lattice-like patterns, 
which phase-lock and travel at constant velocity along 
a common direction. This traveling wave is associated 
with a rectangular lattice, rather than the hexagonal one 
associated with the stationary patterns (Fig. 2g). The se¬ 
lection of a simple traveling wave is favored by periodic 
boundary conditions and small system sizes (SI, Video 
2 ); when choosing large systems (SI, Video 3), or im¬ 
posing no-flux boundary conditions (which might better 
represent experiments) the pattern persists as a flowing 
state, continuously forming clusters of well defined size 
and amplitude (SI, Video 4). 

To further understand the mechanisms of chemorepul¬ 
sive pattern formation, we have performed a linear sta¬ 
bility analysis of our model. Assuming that p relaxes 
fast compared to p and c and F 2 ^ 1 and Dp ^ 1 we 


find p = (s/F)Vc (see SI) and obtain, after rescaling 
variables p' = p/F; d = c/F; n' = /^F, a generalized KS 
(gKS) model [ 12 , 14, 19] that accounts both for chemore- 
pulsion (s = — 1 ), and anisotropic chemical production 
{n i- 0 ): 

p = —sSJ • (pVc) + V^p 

c = skV • (pVc) -f PcV^c p — c. ( 2 ) 

For 5 = 1 , linear stability analysis (see SI) reveals the 
standard chemoattractive KS instability when po > 1 
(or in physical units, when po > [14, 19]. For 

chemorepulsion (5 = —1) we find a qualitatively different 
picture, in which instability arises when npQ > Pc, 01 *, 
in physical units, when ka/3po > ^Dc. In contrast to the 
attractive case, the instability is now determined by the 
anisotropic reaction term; it disappears for isotropic ac¬ 
tive colloids, ka 0. Hence we refer to this as the “Janus 
instability”. In marked contrast to the classical KS case 
for chemoattraction, this is a short wavelength insta¬ 
bility, arising only for q > = ^(1 + po)/(^Po - ^c) 

(Fig. 1 in SI). For large po, or large /^, the corresponding 
length scale is /* ^ (/* ^ (Dpka)/(koVo) in phys¬ 

ical units), and diverges when the isotropic production 
rate vanishes. Remarkably, the steady-state cluster do¬ 
main area grows as oc vq (assuming Dp oc 5]), 

a prediction which turns out to remain valid for Eqs. (1) 
even with finite values of Dp oc v^. This may explain 
recent observations of self-limiting active clusters, whose 
particle number increases linearly with vq [5, 15]. 

Why does the Janus instability lead to arrested, not 
full, phase separation? We find that two effects limit clus¬ 
ter growth. First, each cluster creates a shell of chemore- 
pclient which drives away colloids passing nearby, hin¬ 
dering their arrival. Second, as the cluster increases in 
size, so does the quantity of chemorepellent created at its 
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FIG. 3. As Fig. 2, for s = — 1, po = 8. The snapshots show the amorphous transient dynamics (a), travelling wave patterns (b,c) 
and continously moving states (d) for chemorepulsive colloids resulting from the delay-induced instability. Other parameters 
as in Fig. 2 for (a-c) and Vc — 0.2; = 1; = 0.5; F = 1.0; F 2 = 6.0 and Xu ~ 1.83; tn = 1 for (d). 


core via the isotropic production term: once too large, 
this disintegrates the cluster. 

Note that the gKS model predicts a short wavelength 
divergence of the growth rate, Re[A(q)] as |q| ^ oc, 
suggesting the growth of point-like clusters. We show 
in SI that rotational dynamics destroys this high q- 
divergence. Thus rotation provides an essential ingredi¬ 
ent for chemorepulsive pattern formation, not fully cap¬ 
tured by (2). 

A general, numerical evaluation of the dispersion rela¬ 
tions A(g) 1 ^ 2,3 for the three distinct modes that emerge 
from the full model (1) allows us to plot phase dia¬ 
grams on the plane; a typical example is shown 

in Fig. 4. In general, such phase diagrams show three 
different regimes: (i) one in which the uniform state is 
stable (black in Fig. 4), (ii) one in which the growth rate 
of the instability is real and positive, which corresponds 
to the Janus instability (orange to white); and (iii) one 
in which there is a non-zero imaginary part in the growth 
rate (purple with green stars). Modulo small corrections 
from finite-size effects, the length-scales set by the wave- 
vectors of maximal growth rate closely match those ob¬ 
served in the simulations of Figs. 2 and Fig. 3. 

Remarkably, the phase diagram shows that the oscil¬ 
latory instability can develop even for isotropic chemical 
production, /^ = 0, at large colloidal densities (red ellipse 
in Fig. 4). Accordingly we must have a second instability, 
distinct from the Janus mechanism which requires finite 
K.. (The cyan line in Fig. 4 does not reach k = 0 ai finite 
Pq.) This can be traced to the finite relaxation rate of 
p which we neglected when deriving the gKS model (2). 
Because the oscillatory unstable mode requires in effect 
a delayed reorientation of p, we name this the “delay- 
induced instability”. Intuitively, it can be understood 
as follows (see the simplified, ID cartoon in Fig. Ic). 
Chemorepulsive colloids move towards the minima of an 
initial fluctuation in chemical density (upper panel, solid 
line), accumulate there, and produce chemicals opposing 
the original fluctuation. Due to a finite response time 
this production does not stop at uniform density (mid¬ 
dle panel) but overshoots, leading to reversal and possi¬ 
ble amplification of the initial fluctuation (dashed blue 
line). This cycle repeats (lower panel) and represents a 
delay-induced feedback loop: an initial fluctuation of the 
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FIG. 4. Phase diagram showing instabilities of the uniform 
state for Eqs. (2) and = F = 1. Black denotes 

stability; colors show the wavelength of the fastest growing 
unstable mode and green stars show where this mode is os¬ 
cillatory. The red ellipse highlights instability for = 0. The 
cyan line shows the onset of instability as predicted by the 
gKS model (2). Gircles mark the parameter sets of Fig. 2e-h 
and Fig. 3a-c. Inset: Dispersion relation for the encircled pa¬ 
rameters (a third solutions, not shown, is real and negative). 
Upper arrow heads: fastest growing mode. 


chemical density field triggers another fluctuation of the 
same field but with opposite sign. 

To understand the delay-induced instability quantita¬ 
tively, we develop a minimal model in the SI, showing 
that consecutive fluctuations can amplify if po > 1 lead¬ 
ing to an oscillatory instability. Deep in the pattern form¬ 
ing regime (po ^ 1) we And that the wavelength of the 
fastest growing mode scales as p^ , or in physical units 
as /* oc \/\/3\vQkopo/k^ predicting (in accordance to our 
numerical simulations) that the cluster area grows lin¬ 
early with the self-propulsion velocity. 

The transient or permanent “blinking” of clusters which 
we observed in Fig. 2e appears close to the transition line 
between the Janus instability and the delay-induced in¬ 
stability (regimes (ii) and (iii) in Fig. 1 of the SI). Here, 
stationary and oscillatory modes of different wavelength 
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grow with similar rate out of the uniform state and lead 
to an effective particle motion on top of a stationary den¬ 
sity profile which causes the blinking. (See SI for details). 

In conclusion, our two chemorepulsive instabilities cre¬ 
ate a robust new route to pattern formation. Both insta¬ 
bilities lead to clusters of self-limiting area which grows 
linearly with propulsion speed vq. This agrees with recent 
experimental observations [5, 6, 15] and may shed light on 
the still mysterious mechanism underlying their appear¬ 
ance. (Competing explanations based on the chemoat- 
tractive KS instability either predict macrophase separa¬ 
tion or clusters shrinking with increasing uq [16].) More 
generally, these chemorepulsive instabilities might inform 
design principles for creating active colloids that can self- 
assemble into spatiotemporal patterns with desired prop¬ 


erties. Finally, or key finding that chemorepulsion can 
generate instability of uniform states might also be im¬ 
portant for biophysics, where the chemoattractive KS in¬ 
stability has long been invoked to explain patterns of 
microorganisms [17, 18, 20]. In growing biofilms, for ex¬ 
ample, the interaction of bacteria with self-secreted poly¬ 
mer [24] might be interpreted as chemorepulsion. 
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A. Nondimensionalisation 

To understand the parameter regimes leading to the 
stationary and moving patterns observed in Figs. 2 and 
3 of the main text, we perform a linear stability analy¬ 
sis of our model. We first reduce the parameter space 
by rewriting (1) of the main text in dimensionless quan¬ 
tities (denoted with tildes), i = kgt] x = x^Jkd/Dp 
and p = pko\l3\vo/{klDp)-, p = fo/(yfcd-Dp)p; c = 
vq\P\/{ kdDp). This leads to (now omitting tildes) 

p = -v ■ {pp) + vy 

p = -Fp + VpV^p + sVc - r2|p|^p 

C = PcV^C + p — c + kV • (pp). (1) 

Here F = ^/kd; Vp = Dp/Dp-, = Dc/Dp; k = 
kakd/i^oVo)] F 2 = 72 -Dp/r'o are dimensionless control pa- 
rameters and s = sign(/d) distinguishes positive (5 = 1 ) 
and negative (5 = — 1) chemotaxis. The overall average 
colloid density, po, is a conserved variable and sets the 
initial reference state. Note that the patterning scenarios 
we report in the main text do not change qualitatively 
if this cubic saturation term is omitted and a nonlin¬ 
ear chemotactic coupling is instead used, e.g. 5Vc/|Vc|, 
which prevents unlimited growth of Vc. We favor the cu¬ 
bic saturation term because then r 2 does not enter the 
linear stability analysis of the uniform state. 

As usual, the stability of the uniform state is deter¬ 
mined by the temporal evolution of linearized perturba¬ 
tions around it of wave vector q quantified by their growth 
rates or dispersion relation, A(q). Linearising (p,p, c) 
around the uniform solution (p,p, c) = (po 5 0,po) of (1) 
in one spatial dimension yields: 

/ Sp\ ( Sp\ (dl -Pod^ 0 \ 

{Sp\=M\dp\; M = 0 -F + Vpdl sd, 

\Sc/ ycj \l Kpod^ Vcdl-l) 

( 2 ) 

Inserting a plane wave Ansatz for all three fields (or em¬ 
ploying a Fourier transformation) leads to the following 
characteristic polynomial for A(z), with z \=\q\ 

P(A) := det (M - XI) = {z‘^ - X){VcZ^ - A - 1) 

X {Vpz‘^ — A — T) -h spoz‘^ (^—kX -h z‘^K -h 1 ) (3) 
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As P{X) is a cubic polynomial, we can solve P{X) = 0 
explicitly, leading to closed expressions for the dispersion 
relation A( 2 :). However, the general result for A( 2 :) is a 
complicated expression, whose physical meaning is not 
immediate. 

Thus, we first evaluate the dispersion relations A(2:) 1 ^ 2,3 
numerically which allows us to plot instability or phase 
diagrams (We show and discuss a typical example in the 
main text; see Fig. 4). These diagrams typically show 
three different regimes representing (i) stability of the 
uniform state, (ii) growing stationary modes leading to 
stationary patterns and (iii) growing oscillatory modes 
leading to moving patterns. In the following sections we 
simplify our model to derive simple instability criteria 
allowing for a direct physical interpretation. 


B. Keller-Segel Model and Janus Instability 

To develop a quantitative understanding of the mecha¬ 
nism driving the growth of stationary clusters in Fig. 2e- 
h of the main text (and regime ii in Fig. 4) we now derive 
a reduced model of (1). We first assume that the active 
particle orientation, p, relaxes faster than the conserved 
colloidal density and the chemical density, whose spatial 
average approaches po for long times (c)(t ^ 00 ) ^ po, 
see (1). Hence, we set p ^ 0 and eliminate p adiabat- 
ically. Assuming further weak translational diffusion of 
p, Vp 1, and using that |ppp ~ 0 close to the uni¬ 
form state, we obtain p ^ (5/r)Vc. In this adiabatic 
regime the dynamics constitutes a generalized version of 
the Keller-Segel (GKS) model [1-3], as used to describe 
chemotactic interactions in microbial communities, ex¬ 
tending it to account both for chemorepulsion (when 
s = —I), and anisotropic chemical production (when 

/ 0). In terms of rescaled variables p' = p/T; c' = c/T, 
we can write (now omitting primes). 


p = -sV • (pVc) + 

c = sk\/ • (pVc) + VcV^c X p — c. (4) 


We now perform a linear stability analysis of the GKS 
model. Linearising (4) around the uniform solution 
(p,c) = (po,Po) yields 


(5p\ ^ -spoV^ A fSp\ 


( 5 ) 
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FIG. 1. Dispersion relation for the generalized Keller-Segel 
model for T>c — — 1] s — a — zbl. 


Inserting a plane wave Ansatz and using q := |q| leads 
to the dispersion relation 

2A±(g) =-1 - g^[l + X>c + SK/Oo] ± (6) 

^1 + 2^2 + (2 + k)pqs - 1] + _|_ i^p^g _ i]2 

In the following we discuss the linear stability of the GKS 
separately for s = 1 and s = — 1 by direct analysis of ( 6 ). 
An alternative and somewhat quicker way to derive most 
(but not all) of the following results is possible using the 
Routh-Hurwitz stability criterion [4, 5]. 

Chemoattraction: For chemoattraction or colloids 
turning towards high chemical concentration (s = 1) the 
discriminant in (6) is always positive {Vc^po^K, > 0). 
Thus X{q) is real at all q and we obtain the band of lin¬ 
early unstable modes by solving A(g) = 0. This yields 
q = ±\/(po - i)/(k/)o + ^?c) as the boundaries of the 
band of linearly unstable modes and additionally q = {) 
which is a double zero (We note that A(g ^ Too) ^ 
oc). Hence, linearly unstable modes exist for po > 1- 
In physical units, this instability criterion translates as 
ko/3vopo > Dpkdj; when this inequality is satisfied we ob¬ 
serve phase separation in the GKS model (Eqs. 4). The 
wave number of fastest growing mode follows from the 
conditions dqX{q) = 0 and dqX{q) < 0 as 

_ (1 + ^) y PoKpo + b-l)- b{2po + 6-1) 

9max-± 6(6-1)2 

where b = Kpo + Pc- This wave number determines the 
length scale of the droplets initially growing out of the 
uniform state before coarsening (similar to Fig. 2a, main 
text). 

Chemorepulsion: To determine the instability crite¬ 
rion for 5 = —1 we distinguish two cases. If the discrim¬ 
inant in ( 6 ) is positive, X±{q) is real and the instability 
region follows from X±{q) > 0. In contrast, if the dis¬ 
criminant is negative, the square root in ( 6 ) does not 
contribute to the real part of A±(g) and the instability 
region is determined b —l — q^[l-\-'Dc — i^po] > 0 . In this 
second case, Im(A) 7 ^ 0 and the instability is oscillatory. 

The discriminant vanishes at 

= [a + po A 2^poa] ; a = 1 — Pc + Po(^ + 1) (8) 

and always has a minimum between g+ and g_ (except 
at Pc = 1 + /^Po)- Thus, the discriminant is negative for 
q^e[ql,ql]. 


We can now write the band of linearly unstable modes 
as 

{q\X{q)>0; q^^[ql,ql]}U 

{q\-l-q^[l + V,-Kpo]>0-, y (9) 

and rewrite it (after a lengthy but straightforward case 
by case analysis) in the following form: 

• Regime 1: k. < Pc/po- uniform state stable 

• Regime 2; Vc/po < fc < (1 + V^jpQ + l/p^; 

stationary wave instability with instability band 

> (1 + po)/(/i:po - Pc) 

• Regime 3: k > (l + Pc)/po + l/Po- oscillatory 
wave instability band {q\q^ G ([/^po — (1 +Pc)]~^ 5 ^-)} 
coexisting with a stationary instability band {q\q^ > g?.}. 

Thus, chemorepulsion leads to an instability of the 
uniform state if /^ > Vdpo (in physical units this trans¬ 
lates to: kaPpo > jDc). This instability vanishes for 
K = 0; hence it requires anisotropic chemical production 
by the active colloids and we call it the ’Janus insta¬ 
bility’ (see main text). More specifically, the interplay 
of anisotropic chemical production and chemorepulsion 
generates an effectively negative chemical diffusivity 
(—^cpoV^c term in (4)). If it dominates over the standard 
passive diffusion (PcV^c), this initiates a spinodal insta¬ 
bility. As summarized in the cartoon in Fig. 1 b (main 
text), the colloids, in turn, cross-diffuse towards minima 
in c as a result of self-advection and chemorepulsion 
(the poV^c term in (4)). This cross-diffusion favours 
a colloidal density profile inverse to the chemical field 
(Fig. 2f,g in the main text). 

In contrast to the case of chemoattraction, for 5 = — 1 
the GKS model predicts a fastest growing mode occur¬ 
ring at zero wavelength (A(|q| ^ oc) ^ oc); suggesting 
the growth of very small clusters (compare Fig. 1). In 
the full model ( 1 ) however, rotational dynamics prevents 
this divergence and thus provides an essential feature 
for pattern formation based on the Keller-Segel insta¬ 
bility. This can be seen by performing an asymptotic 
expansion of (3). Rewriting this equation in the form 
P(A) = co( 2 :)-hci( 2 :)A-hC 2 ( 2 :)A^+C 3 ( 2 :)A^ and truncating 
all powers of z but the highest one, separately for each 
coefficient Ci{z) yields 

— A^ J- -f- Pc + 'Dp)X^ 

+ z^^—lXp + Pc(l T T — 0 ( 10 ) 

Solving this equation leads to the following three disper¬ 
sion relations for the three orthogonal modes which are 
linear combinations of the fields p, p, c (these fields form 
a basis of the solution space to ( 2 )) 

Ai, 2 , 3 (g -)> oo) = {-q^,-Vcq^, -Vpq^) ( 11 ) 

Accordingly, rotational dynamics suppresses the Janus 
instability at small wavelength and serves as an essential 
ingredient of pattern formation in our model. We note 
that the Janus instability is robust to sufficiently short- 
ranged repulsions. These interactions provide an alter¬ 
native to rotational dynamics for preventing the small 
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wavelength collapse. 

The orientational order we observed in Fig. 2f can be 
directly understood from ( 1 ) as a consequence of col¬ 
loidal particle number conservation: ( 1 ) imply that any 
current-free stationary state (p = 0, jp = pp — Vp = 0) 
yields p = Vp/p, meaning that every local maximum 
(minimum) of p corresponds to an aster (anti-aster) in 
the p-field. 

Besides unstable stationary modes the above linear sta¬ 
bility analysis of the GKS model predicts also unstable 
oscillatory modes (case 3). It turns out, that these modes 
do not properly explain the moving patterns we observed 
in Fig. 3 of the main text. In particular, Fig. 4 of the main 
text shows that oscillatory instabilities may even occur 
for 5 = — 1 ,/^ = 0 , which cannot be explained based on 
the GKS model. Instead, the moving patterns we ob¬ 
served rely on a delayed response of the active colloids 
to changes in the chemical density. To see this quantita¬ 
tively, we derive a minimal model for this ’delay-induced’ 
instability in the next section. 


C. Delay-Induced Instability 

While the GKS model allowed us to understand the 
Janus instability mechanism leading to stationary pat¬ 
terns, we found above that a distinct instability mech¬ 
anism must be responsible for the moving patterns we 
observed in Fig. 3 of the main text. Specifically, in cer¬ 
tain parameter regimes, a delayed response of the director 
field to changes in the chemical density leads to a qual¬ 
itative failure of the ’adiabatic’ p ^ 0 approximation 
used to derive the GKS model. This delay can cause an 
instability by itself as we discuss qualitatively in Fig. 1 
of the main text. Quantitatively, this can be seen most 
clearly, when switching off rotational relaxation by set¬ 
ting F = 0 (infinite delay). We also set = 0 in ( 1 ) 
and for simplicity Vp = Vc = 0 which yields the fol¬ 
lowing minimal model which is, our simulations suggest, 
sufficient to generate delay-induced patterns of similar 
length scales to the full model: 

p =-V-(pp) + V^p; p =-Vc+r 2 |ppp; c = p-c ( 12 ) 

A linear analysis of ( 12 ) shows that, when /3 < 0 , an in¬ 
stability arises for po > 1 ; in physical units, this reads 
^ol/^l'^oPo > Dpk‘^. Accordingly, fast self-propulsion, 
strong chemotactic coupling, slow chemical decay and 
slow colloidal diffusion all promote the delay-induced in¬ 
stability. This route to pattern formation requires self¬ 
propulsion and works for purely linear reaction dynam¬ 
ics, in crucial distinction to classical reaction-diffusion 
patterning in biological and chemical systems [ 6 , 7]. 

In the minimal model of (12) with infinite delay time 
(F = 0 ) the delay-induced instability first emerges at 
long wavelengths while the fastest growing mode moves 
to shorter wavelength as po increases. We show now, 
that reinstating a finite delay time via nonzero values of 


the angular relaxation rate F generally shifts the insta¬ 
bility to short wavelengths, analogously to the situation 
described for the stationary Janus instability. To this 
end, we reconsider the dispersion relation of the com¬ 
plete model (1) and explore its long wavelength regime. 
The general behaviour of X{q) at small wavenumbers is 
nontrivial. For F ^ 1 (weak damping of long wave¬ 
length modes) however, the fastest growing mode occurs 
at small wavenumbers g ^ l,g 7 ^ 0 (which can be veri¬ 
fied numerically) and we may expand the general result 
for A(g) (the zeros of P(A) given by (3)) first in F and 
then in q. In the special case k = 0, s = —1 (neglecting 
mixed terms of order g^F and higher) this yields 

Re[Xiq)] = -r/2 + q\po-l- Vp)/2 

-q^po{2po-l + 2Vc-Vp)/2 (13) 

for the real part of the the fastest growing mode in (2). 
Hence, the linear instability criterion reads 

(po - 1 - Vp) > 4por[2(po + Pc) - (1 + Vp)] > 0 (14) 

This shows that the delay-induced instability mechanism 
(explained in the main text) is robust with respect to fi¬ 
nite diffusivity of the chemical and director field. Note 
that the instability criterion reduces to requiring that 
Po > 1 when F = Vc = Vp = 0 - this is the exact insta¬ 
bility criterion of the minimal model we just discussed. 

Eq. (13) predicts the fastest growing mode as: 

= Po - 1 - Pp . . 

2po[2po + 2V,-Vp-l] 

Far from the onset of instability (po ^ 1) the wavelength 
of the fastest growing mode scales as pj^^ (assuming 
Dp (X Vq according to ref. Theurkauff et al. in main 
text), or in physical units as /* oc ^|/3|r’o^oPo/^|- The 
selection of a finite wavelength for the patterns requires, 
therefore, a finite chemical decay rate. In particular, the 
length scale of the fastest growing mode scales as I oc yTo 
- therefore the domain area of the corresponding two- 
dimensional clusters scales as oc vq (as in experiments 
[ 2 , 8 ]). 

The wavelength of the fastest growing mode can be 
computed more generally at the onset of instability 
(where Re[A(gniax)] = 0 and we can determine this wave¬ 
length from the condition Re[co( 2 :)] = 0 in (B)) and far 
away from the onset (where A cq and we can neglect 
co{z) in (B)). In general, the scaling for the length scale 
of the clusters, /, varies from I (x Vq (for vq 0 ) to 
I oc Vq (for Vq oo), which means that clusters based 
on chemorepulsion generically grow with increasing self¬ 
propulsion velocity. 

A low-order approximation of the imaginary part of the 
dispersion relation in F and g (neglecting gF and higher 
order terms) leads to ImA(g) = This shows that 

the frequency of the oscillation of long wavelength mov¬ 
ing patterns is given (at lowest order) by cj = In 

physical units, this translates to cj ^ y^(l/6lvokopo)/Dp. 
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FIG. 2. Cartoon of the physical mechanism leading to ampli¬ 
tude blinking of clusters. 


Hence, the frequency of the delay-induced oscillatory in¬ 
stability and the velocity of the travelling wave pattern 
(group velocity associated with the fastest growing mode) 
is, at leading order, independent of the self-propulsion ve¬ 
locity (assuming again Dp oc and the wave number of 
the fastest growing mode g^ax 1/y/^). 

D. Amplitude Blinking 

In the main text, we briefly discussed that amplitude 
blinking of clusters (Videos 1,4 and 5) emerges as a con¬ 
sequence of a simultaneous growth of stationary waves 
with a certain length scale ‘I'Kjqs and oscillating waves 
with another length scale 27r/go and frequency cj. Here, 
we discuss this in more detail and complement the dis¬ 
cussion taking a physical viewpoint. 

Amplitude blinking appears close to the transition be¬ 
tween the Janus instability (regime (ii) in Fig. 4 (main 
text)) and the delay-induced instability (regime (iii)), 
where we have a discontinuous jump of the unstable¬ 
mode wavelength (color change in Fig. 4). The un¬ 
derlying dispersion relation (Fig. 4 (main text) inset) 
shows two distinct, coexisting regions of instability, an 
oscillatory branch at lower q and non-oscillatory growth 
at higher q. Both modes grow simultaneously and the 
fastest growing mode prevails, determining the wave¬ 
length of the emerging pattern. Near the transition, 
where the oscillatory perturbations grow only slightly 


slower than the fastest growing stationary mode, we ob¬ 
serve long-lived oscillations, visible as blinking in col¬ 
loidal density patterns, on top of the Janus instability 
leading ultimately to stable colloidal clusters. The fact 
that in some cases the blinking goes on indefinitely sug¬ 
gests that it can be stabilised by nonlinear terms in the 
equations of motion, not captured in our linear stability 
analysis. 

In the linear regime, the superposition of the fastest 
growing modes of stationary and oscillatory plane wave 
fluctuations creates a chemical density field c(x,t) oc 
[sin(gsx) + cos(cc;t) sin(go^)]- Fig. 2 shows such a field 
at three different times (black lines; for = ^s/2). The 
local minima of this field correspond to the centres of 
clusters (asters) as created by the Janus instability. 

(i) In the upper panel (showing c(x,t = 0)), chemore- 
pulsive colloids move from clusters with comparatively 
high chemical density to adjacent ones with lower chem¬ 
ical density. As we see in the figure, this requires them 
to overcome (small) local density maxima of the chem¬ 
ical density field which is assisted by colloidal diffusion 
(favouring a uniform distribution of colloids), (ii) In the 
middle panel (c(x, 7r/2)), the minima of the chemical den¬ 
sity field are equally deep. The corresponding colloidal 
clusters have different amplitudes (contain different num¬ 
bers of colloids as a result of (i)) but the chemical density 
minima are deep enough to prevent a significant diffusive 
exchange of colloids, (iii) As time evolves the chemical 
density increases at positions of the larger clusters (see 
lower panel, c(x,7r)) which drives colloids again to adja¬ 
cent clusters. This cycle repeats. Effectively, we observe 
a net particle exchange between different clusters result¬ 
ing in the amplitude blinking we observed in Videos 1,4 
and 5. 


E. Videos 

Videos 1-3: Parameters as in Figs. 2e-h, 3a-c and 3d 
(main text) respectively. 

Video 4: Parameters as in Fig. 3a-c but with no- 
fiux (Neumann) boundary conditions instead of periodic 
boundary conditions. 

Video 5: Parameters as in Fig. 2e-h but with Vp = Vc = 
1.104 (This corresponds to the transition line between 
stationary and moving pattern where the growth rates of 
the fastest growing stationary and oscillating plane wave 
modes are identical.) 


[1] Keller, E. F V Segel, L. A. J. Theor. Biol 30, 225 (1971). 

[2] Theurkauff, I, Cottin-Bizonne, C, Palacci, J, Ybert, C, V 
Bocquet, L. Phys. Rev. Lett. 108, 268303 (2012). 

[3] Meyer, M, Schimansky-Geier, L, V Romanczuk, P. Phys. 
Rev. E89, 022711 (2014). 


[4] Hurwitz, A. Math. Ann. 46, 273 (1895). 

[5] Willems, J.L. Stability Theory of Dynamieal Systems (Wi¬ 
ley Interscience Division, New York) (1970). 

[6] Murray, J. Mathematieal Biology. II: Spatial Models and 
Biomedieal Applieations. (Springer-Verlag Berlin), 3rd edi- 












5 


tion (2003). 

[7] Turing, A. M. Phil. Trans. R. Soc. A 237, 37 (1952). 

[8] Buttinoni, I, Bialke, J, Kiimmel, F, Lowen, H, Bechinger, 
C, & Speck, T. Phys. Rev. Lett. 110, 238301 (2013). 



